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In this letter we present a new form of the well known Chaplygin gas model by introducing 
inhomogeneity in the EOS. This model explains u — ~1 crossing. Also we have given a graphical 
representation of the model using {r, s} parameters. We have also considered an interaction of this 
model with the scalar field by introducing a phenomenological coupling function and have shown 
that the potential decays with time. 

' PACS numbers: 

^ ' Recent observations reveals [1, 2] that the present Universe is subjected to an accelerated expansion, 

which can be explained in terms of some new type of matter which violates the strong energy condition 
p + 3p < 0. This type of matter is known as dark energy [3-6], which has the cosmological constant to 
be a strong candidate. However many models have been proposed to play the role of the dark energy. 
Quintessence [7] or the scalar field being one of the most favoured model because of its decaying potential 
term dominating the kinetic term so as to generate enough pressure to drive acceleration. Also one can 
Q^i try Chaplygin gas model [8] with equation of state (EOS), p = —B/p, as it generates negative pressure, 

' where p and p are respectively the pressure and energy density and i? is a positive constant. Subsequently 

this fluid has been modified to p = —B/p" with < a < 1. and 

p = Ap with < a < 1, A, B are positive constants. (1) 

> ■ P" 
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' as generalized Chaplygin gas [9, 10] and modified Chaplygin gas [11, 12] respectively. Modified Chaplygin 

, gas can explain the evolution of the Universe from radiation era to ACDM model. Later inhomogeneity 

^SJ ' has been introduced in the above EOS (1) by considering B to be a function of the scale factor a{t) [13, 

14]. This assumption is reasonable since B{a) is related to the scalar potential if we take the Chaplygin 
Born-Infeld scalar field [15]. 

Interaction models where the dark energy weakly interacts with the dark matter have also been studied 
to explain the evolution of the Universe. This models describe an energy flow between the components. 
To obtain a suitable evolution of the Universe the decay rate should be proportional to the present 
' value of the Hubble parameter for good fit to the expansion history of the Universe as determined by 

^ , the Supernovae and CMB data. A variety of interacting dark energy models have been proposed and 

studied for this purpose [16-19]. 

In this letter we study a new model by considering both A and B in the EOS (1) to be a function 
of the scale factor a{t) and thus introducing inhomogeneity in the EOS (1). We solve the EOS to get 
the energy density and show that the we can explain the evolution of the Universe suitably by choosing 
different values of the parameters. We then consider an interaction between the fluid and the scalar field 
by introducing a phenomenological interaction term which describes the energy flow between them, thus 
showing the effect of interaction in the evolution of the Universe. This kind of interaction term has been 
studied in ref. [20]. 

The metric of a spatially flat homogeneous and isotropic universe in FRW model is 

ds^ = dt^ - a^{t) [dr"^ + r^{de^ + sin^edcj)'^)] (2) 
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where a{t) is the scale factor. 
The Einstein field equations are 



and 



2 = 3^ 



1 



(3) 



(4) 



where p and p are energy density and isotropic pressure respectively (choosing SttG = c = 1). 
The energy conservation equation is 



p + 3-{p + p) = 



(5) 



Now, we extend the modified Chaplygin gas with equation of state (1) such that A and B are positive 
function of the cosmological scale factor 'a' (i.e., A = A{a),B = B{a)). Then equation (3) reduces to, 



p = A{a)p 



B{a) 



with < a < 1 



(6) 



As we can see this is an inhomogeneous EOS [21| where the pressure is a function of the energy density 

p and the scale factor a{t). Also if p = (^^^^ ^^") this model reduces to dust model, pressure being 
zero. 

Now, assume A{a) and B{a) to be of the form 

A{a) = Aoa-" (7) 

and 

B{a) = Boa-"" (8) 

where Aq, Bq, n and m are positive constants. If n = m = 0, we get back the modified Chaplygin gas 
[12] and if n = 0, we get back variable modified Chaplygin gas model. Using equations (5), (6), (7) and 
(8), we get the solution of p as. 



p = a e 



3Aoa- 



Bo f 3Ao{l + a)\'^'^'^^-^"~"' m - 3{1 + a) 3Ao{l + a) 

Co + I I H , a ) 

Aq \ n J n n 



where T{a,x) is the upper incomplete gamma function and Cq is an integration constant 
Now, considering 



P 



for this fluid, we have. 



eff 



AqU " - Boa ''t 



3^0(l + a)a 



Co 



3Ao(1 + q) \— Eq^ C 3^o(l + a)a-" 
n J Aq n' n 



-, -1 



(9) 
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Fig. 1 shows the variation of (jJeff against a{t) for = 1, Bo = 10, a = 1, m = 2, Co = 1 and n = 3(for dotted 

hne), n = 10(for the dark hne). 



where C, = m — 3(1 + a). 

For small values of the scale factor a{t), p is very large and 

B 

p = Ap > Ap 

p" 

where A = Aoa~" is a function of a, so that for small scale factor we have very large pressure and energy 
densities. Therefore initially 

=u}eff = A*a-'' < 1 

P 



where A* is a constant, 



A* = Ao. 



If a = Aq" , the Universe starts from stiff perfect fluid, and if a = 3^o " , the Universe starts from 
radiation era. 



Also for large values of the scale factor 



If 



p = Ap > . 

p" p" 



( = m-3{l + a) <0 



( as we know that upper incomplete Gamma function T{a, x) exists for a < ), the second term dominates 
and hence w — )• —B*a~^, where 



B* = Bo hm e 



3^4o(l + a)a- 



Co + 



3Ao{l + a) \— Bo C 3Ao(l + a)a-" 
n J Aq n' n 



( lima- 



1 and lima^oor(^ 



C 3Ao(l+a)a 



-) — >large value, for ^ < ). This will 



represent dark energy if a > (3^) , ACDM if a = (■^) phantom dark energy 
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if a > ("gr) 3<i+°'"'" . Therefore we can explain the evolution of the Universe till the phantom era 
depending on the various values of the parameters. We have shown a graphical representation of oj^ff in 
fig 1 for different values of the parameters. We can see from fig 1 that Wg// starting from a large values 
decreases with a crosses w = — 1 for some choices of the parameters. 

Since there are various candidates for the dark energy model, we often face with the problem of discrim- 
inating between them, which were solved by introducing statefinder parameters [22]. These statefinder 
diagnostic pair i.e., {r, s} parameters are of the following form: 

r - 1 

and s = — -, p- (10) 



am 3(,-i) 

where H is the Hubble parameter and q (= — p^) is the deceleration parameter. These parameters 
are dimensionless and allow us to characterize the properties of dark energy in a model independent 
manner. The statefinder is dimensionless and is constructed from the scale factor of the Universe and its 
time derivatives only. The parameter r forms the next step in the hierarchy of geometrical cosmological 
parameters after H and q. 



Now, in our case. 



a2 3 

and 



h' = ^, = Ip (11) 



So from equation (10) we get 



a 1 3p 

9 = --E?o = o + o- (12) 



l + -fl + ^)l^---|^. (13) 



2 \ P J dp 2 p da 



(?) 



so that, solving we get. 



! + ?(! + y){Aoa-'' + aBoa-"'x) + hnA^a-" - mBoa-"'x) , s = ^^^^ (14) 
2 2 yy 
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where, y = j = " — B^a "^x and x = p (^+"), p is given by equation (9). 

We have plotted the {r,s} parameters norraahzing the parameters and varying the scale factor a{t). 
We can see that the model starts from radiation era. Then we have a discontinuity at the dust era (for 
radiation era: s > and r > 1; dust era: r > 1 and s — > ±oo; ACDM: r = 1, s = 0; phantom: r < 1). 
The model reaches ACDM at r = 1, ,s = and then crosses ACDM to represent phantom dark energy. 
This model represents the phantom dark energy, whereas. Modified Chaplygin Gas can explain the 
evolution of the Universe from radiation to ACDM and Variable Modified Chaplygin gas describes the 
evolution of the Universe from radiation to quiessence model. 

Now we consider model of interaction between scalar field and the new variable modified Chaplygin 
Gas model, through a phenomenological interaction term. Keeping into consideration the fact that the 
Supernovac and CMB data determines that decay rate should be proportional to the present value of 
the Hubble parameter. This interaction term describes the energy flow between the two fluids. We have 
considered a scalar field to couple with the New variable modified Chaplygin gas given by EOS (6), (7) 
and (8). 

Therefore now the conservation equation becomes 

Ptot + ^^(Ptot + Ptot) = (15) 
so that the equations of motion of the the new fluid and scalar fleld read, 

p + 3H{p + p) = -mp5 (16) 

and 

p^ + m{p^+p^)=2,Hp5 (17) 

( (5 is a constant ). 

Where the total energy density and pressure of the universe are given by, 

Ptot = p + p<t, (18) 

and 

Ptot=P + P4> (19) 

where p and p are the energy density and pressure of the extended modified Chaplygin gas model given 
by equations (6), (7), (8), (9) and p^ and p^ are the energy density and pressure due to the scalar field 
given by, 

p^ = Y + n<^) (20) 



and 



P4> = \-V{<l>) (21) 



where, V{(l)) is the relevant potential for the scalar field 
Thus the field equations become 



= -^Ptot (22) 



and 



^ = -\{ptot + ^Ptot) (23) 
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Fig. 3 shows the variation of V against t for Ao = 1, Bo = 1, a = 5, m = 3, n = 2, Co = 0, S = 0.01. 



Solving the equations we get the solution for p as 



3(l+c«)(l + 5)+n-m 

Bo f 3Ao{l + a) \ " ^^ m-3{l + a){l+6) 3^o(l+a) 

^o + -r 1 z ' z ) 



Ao \ n 



where Co is an integration constant. 

Further substitution in the above equations give, 



(24) 



V{cl,)=3H^ + H+?-j^ (25) 

To get an explicit form of the energy density and the potential corresponding to the scalar field we 
consider a power law expansion of the scale factor a{t) as, 

a = t^ (26) 

so that, for /3 > 1 we get accelerated expansion of the Universe thus satisfying the observational 
constrains. If ^ = 1 or ^ < 1 we get constant and decelerated expansion respectively. 

Using equations (18), (22) and (26), we get, 

3/32 

P^ = ^-P (27) 

where p is given by equation (25) along with (27). Since ptj, is always positive, so we may have at least 
for some range of the values of the free parameters. 



Also the potential takes the form, 



The graphical representation of V against time is shown in figure 3 normalizing the parameters. We see 
that the potential decays with time. 

Here we present a new variable modified Chaplygin gas model which is an unified version of the dark 

matter and the dark energy of the Universe. It behaves like dark matter at the initial stage and later it 
explains the dark energy of the Universe. Unlike the Generalized or Modified Chaplygin gas model, it 
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can explain the evolution of the Universe at phantom era depending on the parameters. Also we have 

calculated the {r, s} parameters corresponding to this model. Normalizing the parameters such that 
m — 3(1 + a) < 0, show the diagrammatical representation of {r,s} for our model (in Fig. 2), varying 
the scale factor. We see that starting from the radiation era it crosses ui = —1 and extends till phantom 
era. Also we can see that the deceleration parameter starting from a positive point becomes negative, 
indicating deceleration initially and acceleration at later times. Again we have considered an interaction 
of this fluid with that of scalar field by introducing a phcnomenological coupling term, so that there 
is a flow of energy between the field and the fluid which decays with time, as in the initial stage the 
fluid behaves more like dark matter and the field that of dark energy, whereas in the later stage both 
explain the dark energy present in the Universe. In Fig. 3, we have shown the nature of the potential by 
considering a power law expansion of the Universe to keep the recent observational support of cosmic 
acceleration, and we see that the potential decays with time. 
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